Topological Expansion for the Cauchy two- Matrix- Model 



M. Bertolat'ffl A. Prats Ferrer^ 

* Department of Mathematics and Statistics, Concordia University 
1455 de Maisonneuve W., Montreal, Quebec, Canada H3G 1M8 
" Centre de recherches mathematiques, Universite de Montreal 

2920 Chemin de la tour, Montreal, Quebec, Canada H3T 1J4- 
Abstract 

Recently a two-matrix-model with a new type of interaction [T] has been introduced and analyzed using bi- 
orthogonal polynomial techniques. Here we present the complete expansion for the formal version of this 
model, following the spirit of i.e. the full expansion for the non mixed resolvent correlators and for the 

free energies. 
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1 Introduction 

The relationship between matrix models and graph counting has attracted attention for a long time. Such 
a connection hinges on the interpretation of Feynmann diagrams in the formal expansion as random lattices 
with or without matter (see for small selection of important works). Very recently the full structure 

of this expansion has been unveiled [HI [21 El E] • The construction relies on the definition of a formal matrix 
model that provides the generating function for counting of fat graphs embedded in Riemann surfaces, and the 
topological expansion is the organization of these by the Euler characteristics of the Riemann surface in which 
they are embedded. The first type of matrix models where this structure was completely revealed was the one 
(Hermitian) matrix model, which counts orientable graphs with no other decoration [10l [3] . 

In multi-matrix models usually only the polynomial interaction in the action is considered. Apart from some 
isolated works, most of the attention has been directed towards the so-called Harish-Chandra-Itzykson-Zuber 
type interaction which appears to be the only generically integrable casjj]. The first multi-matrix model for 
which the full topological expansion was revealed was the two (Hermitian) matrix model, and happened to 
formally share the same structure [H [9] . Similarly, in [11] the same structure was found describing the full 
topological expansion of the chain of (Hermitian) matrices. 

This apparent universal description relies only upon the existence of an algebraic curve describing the leading 
order of the topological expansion. In [12] the same construction was developed for any (non-singular) algebraic 
curve without considering any underlying matrix model. The construction has been successfully applied in 
connection to certain non-algebraic curves and in relation with algebraic-geometric counting problems (see 
[TJirriJIin] as an example). 

4 In some sense it is the only Gaussian interaction possible. 
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Although the general case has already been worked out, the combinatorial interpretation is not clear for an 
arbitrary algebraic curve, except for some isolated particular cases. 

In this work we show how the same topological expansion construction works in a recently introduced two 
(Hcrmitian) matrix model called Cauchy two matrix model pQ. 

In section [5] the Cauchy matrix model is recalled. The eigenvalue representation is presented and then 
written back into the matrix representation in a much more convenient form. 

In section [3] the two master loop equations needed for the construction are presented. The main difference in 
the derivation (contained in App. [Aj is that the domain of integration has a natural boundary since it extends 
only over the positive definite matrices. This yields some poles in the structural master loop equations. 
In section |4] the algebraic curve (first step of the construction) is found as the leading order contribution from 
one of the master loop equations. This will provide the leading order solution for the one point resolvent. 
Everything is expressed in terms of the algebraic curve. 

In section [5] the other master loop equation's topological expansion is analyzed. This leads to a recursion 
equation between the different orders in the topological expansion of 1 and 2 point resolvents. The recursion 
equation is solved and the solution for general fc-point resolvents is presented. 

In section [6] the topological expansion of the free energy (thus that of the partition function itself) is found by 
inverting the loop insertion operator on the one point correlator. All the terms except the order zero and the 
order one can be found in that way. 

In the appendices we have detailed some of the calculations that do not bring great inside into the problem or 
that are fundamentally equivalent to similar calculations in other works. 

In|B]we give the definitions of the objects related to the algebraic curve that we use, and give some relations 
and properties. 

In [C] we compute the order one term in the topological expansion of the free energy that could not be found by 
inverting the loop insertion operator. The computation is relatively standard and it is included for convenience 
of the reader. 



2 The Cauchy Matrix Model 

Consider the following matrix integral 



e -f tr(Fi(Mi)+V 2 (Af 2 )) 

Mi=Ml>o ' ' det (Mi + M 2 y 



Z c = I d.\/id.\/ 2 — - v <• •••• / ' (2.1) 



where the domain of integration (as indicated) is the ensemble of N x N Hcrmitian positive definite matrices. 
We will view the matrices Mj as random variables with the probability measure 



dP(Mi,Mi) := — dMidM 2 r - 77Z —rr^f (2.2) 



1 e"T tr(Vi(Mi) + V 2 (M 2 )) 

— dMidM 2 

Zc det(Mi+M 2 ) 

The potentials Vi parametrize the measure associated to each matrix, and must satisfy a suitable growth con- 
dition at infinity for the integral to be convergent; as typical we will require Vi(x)/ln(x) — ► +oo. 

x — >+oo 

We will also assume that they diverge at the origin of the spectra with at least logarithmic growth. The 
simplest class of such potentials, to which we restrict in this paper, is the class of "polynomials + logarithm" 

dk 

Vi{x) = -t { \hix + ^tfx k , >0,t ( l\>0 , j = 1,2 (2.3) 
fc=i 
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although later we will allow (formal) variations with respect to the infinite number of parameters tl . The log 
term is included to avoid the eigenvalues approaching the hard edge of the spectrum of positive matrices at 
x = 0. The expression T is called the free energy of the matrix model and we call T the total charge. 

Thanks to the Harnad-Orlov formula [TC] and following [T] we can write (|2.1[) in eigenvalue representation 

Z C = [Yldx^yA A 2 ( Xl )A 2 ( yi ) (2.4) 

where Xi (resp. yi) are the (positive) eigenvalues of M\ (resp. M2), and A(a^) = det(a^ _1 ) = 

V / l<iJ<N 

Yii<j( a i ~ a j) i s the Vandermonde determinants. 

Note that we can easily write (|2.4[) back into matrix representation with a different (but equivalent as long 
as we look at non-mixed correlators) interaction, namely 

„-f tr(V 1 (M 1 )+V 2 (M 2 )) 

dM 1 dM 2 : r (2.5) 

Ah>0 det(Mi®I + I®M 2 ) v ; 

where now the determinant is acting on a TV 2 x TV 2 size matrix, and I represents the TV x TV identity matrix. 
Note that the TV exponent in the determinant from (|2.ip has disappeared. 

We will consider here the formal perturbative expansion of the matrix integral around one local minima. 

The goal of the following sections is to analyze the loop equations for this matrix model and solve them in 
the spirit of [2]. For this we have to restrict ourselves to the formal version of the model, i.e. by choosing one 
of the many local extrema of the potentials and perform a perturbation of the action around it. This is in fact 
equivalent to the fixed filling fraction condition that amounts to specify and fix the amount of eigenvalues that 
lie close to a given extrema of the potential. Indeed a local cxtremum for the matrix action can be characterized 
by the number Tvj of eigenvalues of the matrix Mk (k = 1,2) that lie close to a given local extremum j for 
the eigenvalue action. This numbers (after normalizing them) are called filling fractions. All this considerations 
assure the existence of the topological expansion that we are trying to find. Equation (|2.5p is the most convenient 
form of Zq for this type of analysis. 

3 Loop equations and spectral curve 

This matrix model is somehow peculiar because instead of one master loop equation we will need two of them. 
One of them will provide the large TV algebraic curve that underlie all our calculation. The other contains, in its 
large TV expansion, the recurrence equations for the algebro-geometric objects that are naturally associated to 
the algebraic curve and that solve the model. Since the derivation of the loop equations is a somehow standard 
technique, the ones we use in this paper are derived in appendix [XJ In this section we only define the objects 
appearing in these loop equations and write down the equations before analyzing and extracting all information 
we need from them. The interested reader can find more details on their derivation in appendix IA"! 

3.1 Definitions 

Let us first define the following scalar random variables depending on the spectral parameter x, 

T ( 1 
Wi(x) = -tr 



N \x-Al! 



TV V x - Mi 

t 2 f ^-vgm) ( 



TV 2 V a; -Mi lMi®I + l®M 2 
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and of M 2 and the same spectral parameter x 

T ( 1 
W 2 (x) = --tr 



W N \ x + M 2 J y ' 



N 2 V x-Mi \Mx®I + I<g> M; 

The notation tr, (...) stands for the partial trace on the corresponding factor of the tensor-product space. 
The functions Wi(x) are called resolvent and will be considered only as formal expansions in x , i.e. as a 
generating function for traces of powers of Mj. The functions Pi and Qi are meromorphic functions of x with 
a pole only at x = 0. 

Finally we introduce the following x-independent scalar random variables 



h = ^tr(Vl(M 1 )) + ^tx 1 , 2 r 



N v iV " N 2 ' \Mi®I + I®M2 

T T 2 ( 1 

tr(^ , (M 2 )) + - J tr 1 



(3.3) 



v 2V " iV 2 ' \M 1 ®I + I®M 2 
The dependence on Mj of 7^2 will be understood. We introduce the following notation for the average 



/ \ I f e -^tv(V 1 (AI 1 )+V 2 (M 2 )) 



>c J m 4 >o det (Mi (g) I + 1 (g) M 2 ) 

Finally we define the following auxiliary functions for later convenience in writing the loop equations 

2V{(x) - Vi{-x) 



Y 1 (x) = Ufa)-W 1 (x), Ufa) 
Y 2 (x) = Ufa)-W 2 (x), Ufa) 



3 

~V{{x) + 2V 2 '(-x) 



(3.5) 
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They differ from the corresponding random variables only in a deterministic shift that depends explicitly on 
the potentials. Although it may seem redundant we also need the following definitions 

Wfa) := -Wfa) -Wfa) 
Y {x) := -Yfa) - Y 2 {x) = Ufa) - Wfa) 



(3.6) 



3.2 The Master Loop Equations 

The first of the two master loop equations is quadratic in both Y% and Y 2 and has the form 

( (Yfa)) 2 ) + ( (Y 2 (x)f) + (Yfa)Y 2 (x)) = (r(x)) ~l(h+ h) := (£(*)) (3.7) 

where 

R(x) = \(V{(xf + Vfaxf - V{{x)Vi{-x)) - P 1 (x) Pfa) 

= (Ufa)) 2 + (Ufa)) 2 + Ufa)Ufa) - P 1 (x) - P 2 (x) 

is a meromorphic function with poles of order at most 2 only at x = 0. 

The second of the two master loop equations is a cubic equation that both Y\ , Y 2 and Y3 satisfy identically 

(%(*)») - (R(x)Y k (x)) - ^ + '-I) (Wfa)) = (D(x)) , for k=l,2 (3.9) 
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where 

[D(x)) = -U' Q {x)U[{x)U' 2 {x) U{(x) (P 2 (x)) - U' 2 {x) - \ (u[(x) (j 2 ) + ^ (h)) 

+ (Qi(x) + Q 2 (x)) + (W^x^ix)) + (w 1 (x)W 2 {xf) + (Si(:r)) + (s 2 (x)) 



(3.10) 



is a meromorphic function with pole at x = 0. The last two terms are defined in the appendix [XJ The last four 
terms contribute only to the residue at x = 0. 

4 The cubic master loop equation and the spectral curve 

Before going any further let us rewrite equation (|3 . 0|) in a suitable form for latter use. As usual with matrix 
models loop equations we rewrite everything in terms of connected correlators. 

Let A be any scalar random variable obtained as trace of a polynomial of Mi . We define the connected part 
of a two-trace correlator as 

iV 2 // \ / \ / \\ d / dA(£) 



A(0w 1 (O) c = - 2 {(a® Wi ®) - (a( ) (w l( o)) = ^ (a( ) - (^) (4.1) 



with the definition of the vertex operator 

This formula is valid even when A(£) depends on V\. 

Remark 4.1 Note that the last equality reduces to the usual definition whenever A{^) is independent of the 
moduli in V\. When A depends explicitly on V\ we need to add the extra term in the last equality of (|4.1|) to 
compensate for this dependence. 

In particular we are interested in the following relations 



(Y k (0f) = (n(0) 2 + ^ (w k (Ow k (o) c 
(r k (0YA0) = (r k (0) (m)) + ^ (w k (Ow k (o) c ( 4 - 3 ) 



where we have defined 



R(0 = -Pi(0 - p-2(0 - j{h + h) (4.4) 



Similarly we obtain 

( (Y k m 3 ) = (r fe (0) 3 + ^3 (Y k (0) (w k (0W k (0) c ^ (w k (0W k (0W k (0) c ■ (4.5) 
Putting all together we have that ()3.9j) becomes 

Y k (0) 3 (R(0) (Y k (0) = (d(0) + ^ {w k (0 3 



(4.6) 
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4.1 The large N limit of the cubic master loop equation 

Wc now postulate a T 2 /N 2 expansion of the fornjfl 



oo , T s 2h 

<n(o)=»(o-E(4) < h) ® 

h=i ^ ' 

(a(o) = # 0, (fl+£(|) 



(4.7) 




2/1 



and take the leading term of 14.61 



(4.8) 



The large iV limit of the cubic loop equation gives us the spectral curve 



E(x, y)=y 3 - R (0) (x)y - C (0) (x) = 
whose solutions are encoded in the large N limits of ( Yk(x)j 



(4.9) 



y = y k (x) := Urn (Y k (xj) , k = 0, 1, 2 



(4.10) 



The relation Sj^offcO 15 ) = (which is obvious from the definition of Yq(x)) is the reason of the absence of the 
quadratic term in the cubic spectral curve ()4.9|) . 

The meaning of Ansatz (|4. T[) -as usual in these computations- is that it is consistent with the formal 
manipulations of the loop equations, which is our perspective in this work. The expansion for the large N 
expectations in the convergent model typically will not be of the advocated form, in particular if the asymptotic 
spectral curve (|4.9p has a genus other than zero. The fixed filling fraction condition placed us from the beginning 
away from the convergent model (except for the genus zero case also known as one cut case) . 

4.2 The spectral curve 

We now analyze the basic structure of the spectral curve E(x, y) that we denote by S. In terms of the x variable 
S is a three-sheeted covering of C. Each sheet is parametrized by one of the solutions yu{x) of the cubic equation 
(|4.8p . It was shown in [TJ [T7] (in the context of the convergent matrix integral) that the sheets 1 and are glued 
along one or more cuts contained in R+ and sheets 2 and along cuts contained in R_. 

We will think of the three functions yk{%) as one single valued function y{p) on S where p € S. Similarly 
we think of a; as a function x : E — * P 1 so that E(x(j>), yip)) = for all pgE. Let us call \k the sheet defined 
by yk(x), i.e. 



In general we will always view all the functions previously introduced as single-valued on the spectral curve, 
for example, the ones collecting Wk{x) and U' k (x) 



Xk = such that y(p) = yk(x(p))} 



(4.11) 




whenever p G Xk 



whenever p e Xk 



(4.12) 



5 Indeed, by definition the 1/N 2 expansion of / ) is the same as that of Wk(£,) except for the leading term. 
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Since the projection x : £ — > P 1 
is three-sheeted, there are three different 





points pW (i = 1, 2, 3), such that x(jp^) = 
x(p^>) for all i and j. Whenever we write 
superscripts to the points on the curve, 
we will be referring to points on the same 



/ 



■X 




^-projection. There are three different Figurc L Thc choicc of homology basis . The dashed contours run 



points above x — oo for y(p), one in each n xi, the solid contours on X2 and the dotted contours on xo- 
sheet. We will call them oofc, and in a 

neighbourhood of ook we have the behavior y(x) ~ U' k (x(j>)). Other important points are the endpoints of the 
cuts. These are solutions to the equation dx(p) = and are the only points belonging more than one sheet. We 
assume that they are simple zeros of dx(p) which means that these points belong to exactly two sheets. 

We denote the endpoints of the cuts by a- fe ^ (k = 1, 2, i = 0, . . . , 2d~k — 1, with dk < dk) where the superscript 
indicates whether they belong to xi H Xo or X2 H Xo- Assume that |aj) | < \a^\ < ■ ■ ■ < |<4j _il> then tnc cuts 
are the intervals A^ = [oS^ , c4i+i] for i = 0, . . . , dfe — 1. We will call on a generic branch point when there 
is no need to specify to which family it belongs. 

We must also introduce a canonical choice of homology basis for S. Not all the cuts of each species are 
independent since the sum of them all is reducible to a point. There are dk — 1 independent cuts of the 
species^], thus the maximal genus for S is d\ + d~2 — 2. We define the cycles A^ on xi encircling only the cut 
A^ for i = 0, . . . , d\ — 1 and similarly the cycles A^ on X2 encircling only the cut A\ 2 ^ for i = 0, . . . , di — 1, 
and choose as independents all except A^ = — J^'j'Li 1 ■Aj (k = 1, 2). The conjugated cycles go through 
the 0th and the fcth cuts on Xk visiting both Xk and xo- It is easy to verify that those cycles form a canonical 
basis (see Fig. [TJ. 

4.2.1 Moduli of the model. 

The data that defines the model is contained in 

( k) 

• the values of the potential parameters tj ; 

• the total charge T; 

• and the filling fractions defined below in eq. (|4.17[) . 

We denote a generic setting of these parameters as M = {tj , T, tf^}- We will be interested in computing 
variations (derivatives) with respect to these parameters of different functions of them. In particular, for the 
potential parameters, we introduce a generating operator of such variations. From the matrix model point of 
view the so called loop insertion operators are an important object. They define the one point resolvent when 
applied to the free energy, and define the connected n + 1-point resolvent when applied to the connected n-point 
resolvent. In our model we have two of them: one for each matrix 




(4.13) 



Remember that dk is the degree of the polynomial potential V^(x), and that dk < dfe. 
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With these two operators we will define the combined loop insertion operator in £ 



a 



whenever p £ xi 



- — whenever p e \o (4-14) 



dV( P ) 



dVi(x) dV 2 (-x) 
d 

dV 2 {-x) 



whenever p G \ 2 



We remark that from the above definitions it follows by a direct computation that 

d 



av(p) 



T = W(p) 



~ W(p) U (P ' " 3 (x{p) - x(p')) 2 
4.2.2 Fundamental differential on the algebraic curve 

The curve £ can be parametrized equally well by Ai or by the meromorphic functions y(p) and x(p) (defined on 
some union of open sets in C). In fact we will see that the differential y(p)dx(p) on £ contains the information 
in A4. Indeed,we can extract vjp for k > from the behavior of ydx around_| 0012 

Rcs y( P )dx(p) = 2tg - (-l) fc 4 2) 
Res y( P )dx(p) = + (-l) fc 2t< 2) 

p— >00 2 3 

Also the filling fractions can be explicitly written as the A cycles of the fundamental differential y(p)dx(p) 

— i y(p)dx( P ) = ef = (-l) fc J ' ■ (4.17) 

The shifts ?/ fc ) m the total charge T above are due to the residues t_\, k = 1,2 at of the derivatives of the 
potentials and are related to those as follows, 

2t {1) +t {2) 

2^+t« (418) 

(2 ) = + t_i 

' 3 

The sign in (|4. 17)) appears because of the reflection to the negative axis. The filling fractions condition 

X^o 1 e j = (-l) fe (T + rf k ^) f° r both- k = 1,2 is another way of writting Aq — — Y^j=^ Aj and implies a 
relation for jyW and the total charge T 



Res y(p)dx(p) = — 77^' + rf 2 ' 
(i 



Res y( P )dx(p) =T + if L) (4.19) 



J) — >oo± 

Res y(p)dx(p) = - (T+ 77 (2) ) 

p— >002 



4.2.3 Multidifferentials on the algebraic curve 

In the next two sections we are going to compute the full topological expansion for the n-point resolvent 
correlator and the free energy, and in order to do so we need to define some functions and differentials on the 
curve. Consider the n-point resolvent correlator 

W ku „„k n (x(pi),...,x(p n )) = (f[W ki (x(p i ))\ . (4.20) 



7 The behavior around the third infinity point p = ooq is not independent due to the relation 3/1 (x) + yiix) + 1/3 (x) = 0. 
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the construct a function on E for every n 



W n {p u ■ ■ ■ ,Pn) = W kl _ kn (x( P i), ■ • .,x( Pn )) fo^iH = I' ! X n. (421) 



and its topological expansion 

o 

W n ( Pl ,..., Pn ) = J2(^) Wi h \ Pl ,..., Pn ). (4.22) 



rp\ 2h 



h=0 

Finally, take /i-order n-point resolvent correlator and construct the multi-differential 



5jW(pi, ..., Pn ) = Wi h \ Pl , . . . , Pn )Y[dx( Pl ) (4.23) 

i=i 

In particular we have the one-point differentials 

uj [ l l \p) = W W (p)dx(p) (4.24) 

Remark: Note that the notation changes and the subindex no longer indicates which sheet the point lies. Instead it indicates that 
it is a n-differential or n-variables function. 

Another way of writing it using the loop insertion operator is 

n n o 

sjWfa, . . . , Pn ) = n mp 1 ) n whA^ (4 - 25) 

i=l i=l ^P 1 ' 

where is the hth term in the topological expansion of the free energy J- . 
For latter use we define also 



!h)t \ -(h), \ i x x dx(pi)dx(p 2 

(x{pi) - x(p 2 )Y 



(4.26) 



which only modifies o;^ "* (p) and leaves the rest untouched. 
4.3 Variations of the spectral curve with respect to the moduli 

We begin by studying how y(p)dx(p) changes under variation of the moduli. Since we have the constraint 
E(x,y) = 0, x(p),y(p) and p are not independent, we have to keep one of them fixed while computing the 
variations. For convenience all variations are taken at x{p) fixed. We will generically call this variations dQ. 
The variations will be written in terms of standard objects in algebraic geometry that we define in appendix [Bl 
The reader can find there all necessary information. 

4.3.1 Variation with respect to the filling fractions 

Not all filling fractions are independent due to the restrictions pointed out in the previous section. As with 
we consider as independent the for j = 1, . . . , d k — 1. Clearly, since these ones are independent we 

have — fjry = S it jSk,k' for i, j ^ 0. 

Let us compute the variation of y(p)dx(p). We have that 

V y (p)dx(p) = -4y(^(p) - W{p))dx(p) = O {x{p)- 2 ) dx( P ) (4.27) 
de) ' de) ' 
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when x(p) — > oo. Indeed the potential U'(p) does not depend on e's by hypothesis, and W(p) gives us the 
behavior at oo. On the other side we can compute the A cycles of the above differential 

2^£w W) v{p)dx{p) = ^T2^£w yip)Mp) 

All together we have that - ^ fc) y(p)dx(p) is a differential on £ with no pole and with normalizing condition 
(|4.28[) . which determine them completely to be a basis of normalized first type Abelian differentials defined in 
appendix IBl and that can be written in terms of the Bergmann Kernel B(p, q) defined as well in appendix [51 

dn ik) -^y( P )dx(p) = 2mduf ] = I B(.,p) (4.29) 



- a (fe) ; 

4.3.2 Variation with respect to r/ fc ' and the total charge T 

Imagine now that we change the total charge T, the variation of y(p)dx(p) will be 

{O (x(p)- 2 ) dx(p) pe X o 
(7(k + ( x (p)~ 2 )) dx (p) P^Xi (4.30) 

Indeed by hypothesis the potentials do not depend on T, and the behavior at infinity of Wk(x) gives us the 
above behavior. Clearly, this differential has poles only at 001,2 of order 1 with residues +1, —1 respectively. It 
is also easy to see that the A cycles of this differential are zero because the filling fractions do not depend on 
T either. The conclusion is that 

dn T ( P ) = —;y(p)dx(p) = dS 001 ,o 0a (p) = / B(. :P ) (4.31) 
which is a normalized third type Abelian differential (see appendix [Bj . Analogous calculations show that 



d 



d ^(i)b) = J^(T)y(p) dx (p) = dS oo u oo (p) = / B(.,p) 

'000 



d 

d^( 2 )(p) = Q^2jy(p) dx (p) = dSoo ,oo 2 (p) 



X. 1 



(4.32) 



B(.,p) 



X 2 



4.3.3 Variation with respect to the potentials 

Consider now the potentials. The variations with respect to the parameters of the potentials look like 

dfi t ci>(p) = ^y(p)dx(p) = —^(U'(p)+Tx(p)- 1 +0{x(p)- 2 )dx(p) 

2 



jx{p) J 1 dx{p) + O (x(p) 2 ) for p e Xi 
L yx(py- 1 dx(p) + O (x(j>)- 2 ) for p <£ xi 

dn tW (p) ee JLy(p)dx(p) = -4y(C/'(p) - Txip)- 1 + O (x(p)- 2 ) dx(p) 



(4.33) 



r_i v- 1 I ii x (py~ ldx (p) + o {x( P )- 2 ) for P e X 2 

{ ' \ -\jx{py- 1 dx{p)+0(x{p)- 2 ) forp^ X 2 
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This behaviors at oot plus the normalizing conditions j>Ak) dQ, J 
nation of second kind differentials (see appendix [B| 



(AO 



determine completely dfi^ as a combi 



dOf (p) = - Res x(.yB(.,p) 

O 00 1 



- Res x(.yB(.,p) 

J 00 2 ,000 



dOf(p) = (-iy- 1 (|Re Sa; (.) J 'S(.,p)-i Res x(.yB(.,p) 

\ J 002 O 001,000 



(4.34) 



Explicit vs implicit derivatives of ydx. We claim that we can represent ydx 



2 d k 

ydx =2ttz^ ^ ej 



TdS'o 



(p) + ?7 d-^ooi.ooo + »? 



(2) ds 



fe=l j=l 



fe=i j 



(fe) 

i 



(4.35) 



where i a is a collective notation for the whole set of parameters of our model Ai = {tj , T, 6j , r]^}. To see 
this, we provisionally denote by 77 the right-hand-side of (|4.35[k it is seen, from the definition of the various 
terms that 77 — ydx has no poles at 000, 00 1, 003 and hence it is a holomorphic differential. On the other hand, 

(k) 

from the definition of the filling fractions e^- it also appears that it has vanishing „4-cycles. Standard theorems 
then guarantee then that 77 — ydx = 0. Therefore (|4. 35)) yields an alternative representation of the fundamental 
differential ydx. 

On the face of it appears that dt a ydx is the same as the coefficient of the corresponding "time" t a in the r.h.s. 
of (|4.35p . namely, the "explicit" derivative of the r.h.s. of (|4.35p . Note that this is not a trivial statement, since 
the spectral curve and all the differentials appearing as coefficients of the times in (|4.35[) do depend (implicitly) 
on the times as well. 

Note that all variations df2 ta are expressible as an integral operator on the Bcrgmann kernel B(p, q) 

dO t _ = —y(p)dx(p) = I B(p,-)A ta (-) (4.36) 



dt a 



where C to is the integration path associated with the variation with respect to t a and A ta (p) is a meromorphic 
function on the curve. For example, looking at the formulae above, we have C t i — ^C OQ1 — i (Coo + Coo 2 ), an( l 

A.(d(p) = x(j>y , where Coo k is a small cycle around oofe. We thus introduce the integral operators associated 

j 

to the various derivatives 



dt 



(1) 



dt 



(1) 



□ 



□ 



□ 



d 




□ 



dT 



□ 



2 1 

- Res — Res ) x J □ 

3 001 3 000,002 j 

2 1 

— Rcs+- Res 1 (~ x ya 

3 °°2 3 °°0 1°°! 



□ 

OO 1 



(4.37) 



□ 



300 
000 



□ 



CSOi 



□ 
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where the regularized integrals arc defined (if necessary) by subtraction of the singular part in the variable 
x(p) (see [18]) With this notation we have -compactly- 

y(p)dx(p) = Y, t a Ja(B(;p)) (4.38) 
a 

4.4 The two-resolvent correlation function 

We can now compute the action of the loop insertion operators on y(p)dx(p). 

Notice that a part from a constant factor in front depending on which sheet p sits, the sum in the above formula 
produce a second order pole whenever x{p) = x{p'). This fact, together with the vanishing A periods condition, 
indicates that 

d 2 



Mp') 7 ^y(p)Mp) = ~B(q^,p) + \ (B( q ( \p)+B(qW,p)) 

aV^ (x) 3 3 V / , 

(4.40) 



3 (x(q) - x(p)) 2 

where B(q,p) is the so called Bergmann kernel, and q^ £ \k satisfy x{q^ k >) = x(q) = x. In the last line we 
have used the relation 

V B(qW p) = dx ^ dx ^ (4 41) 

We obtain similar results by applying dx(q) Qv ^ x{q}) or -dx(q) ( dVl d x ( q)) + ov 2 (-x( q )) )- AU thrcc operators 

d 

dVjg) 



combined in dx(q) g y, > give 



where d/dV(x) was defined in P~H| and U in (|4~T2")) . 
It is now a simple calculation to find that 

u^\q,p) = dx(q)^-4 0) (p) = dx(q)^- ) (U'(p) - y(p)) = 

-B(q,p) dx (l)Mp) (4.43) 

(x(q) - x(p)) 2 
= -(B(q+,p)+B(q-,p)) 

where q^ are the two gW ^ q. Note that (|4.43[) implies (p',j») = B(p',p). We can also write it as 

c4V,P (fe) ) = -E^°V,P W ) (4.44) 

which will be used latter. 

5 The quadratic master loop equation 

Once we have extracted the algebraic curve from the cubic loop equation, only remains to find the tower 
of recursion relations deriving from the master-loop equations. The recursion relations for the terms of the 
topological expansion are contained in the topological expansion of the quadratic loop equation. 
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5.1 Recursion relations and topological expansion 

Using (|A.26[) we can rewrite equation (|3 . T[) as 



J2(Y k (x) 2 ) = 2(R(x)). (5.1) 



fc=0 

Indeed, introducing the definition (|A.26[) in the LHS of (|5.1|) we get twice the LHS of (|3 . 7[) . Write this in terms 
of connected components 

2 2 

E <^^)> 2 - 2 - ^ E (^) 2 } c ( 5 - 2 ) 

fe=o fc=o c 



and extract the order — jtt of t ne previous equation. We get the equation 

2E»(P W )W 1 W (P W ) = - 2i? ( ' l) (a;) + £ f]T ^"V^f"™^) + Wf^W^ (5.3) 

fc=0 fe=0 Vm=l / 

for h>0. 

Equation (|5.3|) is of the same form as the recursion equations solved in [2] and generalized in [12] to arbitrary 
algebraic curves. The final solution is written in the following way 

and can be generalized to 

wj&fapK) = ~ E 5 (y(p) d -7(P))dx(p) S w iTi(P.Pj) w fc-j+i(P.PK\j) +4+2 X) (p.^Pk)^ (5.5) 

for higher correlators by applying to equation 15.41 The notation j>k stands for an array of variables 

Pit,Pi 2 , . . . ,Pi k where K = {ii, i2, ■ ■ ■ , ifc} is the set of subindices, J stands for subsets of subindices of K, and 
K\J is the complementary subset of J. The derivation of these identities is contained in [2] and we refer ibidem 
for details. 

Remark 5.1 In those recursion relations, the initial data is contained in y(p)dx(p) and oj^ (p,p'). Note that 
the latter is defined by (|4.26[) and (|4.43[) to be ui^\p',p) — B(p' ,p). It is not exactly equal to the 2-point 
resolvent correlator (remember the shift with respect to TD^ )■ On the other side, the other differentials cjjj. d 
coincide with k-point resolvent correlators. 

Remark 5.2 It is also important to notice that by the differentials dS(q) in those formulas are defined through 
integral formulas, e.g. 

dS PtP (q)= f B(q,r) (5.6) 
Jp 

where the whole integration path must remain in a small neighbourhood of the corresponding a, or be homolog- 
ically equivalent. 
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5.2 Variations of ui^ with respect to the moduli 

In [12] the authors proved (see section 5.1 and Theorem 5.1) the following relation 

W a W * ){pK)= l c J&facqjAtM (5.7) 

for differentials satisfying (|5.5p and variations of the form (|4.36l) . This equation is in fact essential in order to 
find the topological expansion of the free energy. The proof is by induction supposing that it is satisfied for 
all w[ up to some (fc, /i)-level. This important formula is telling us the way these correlators change when we 
change our model. This formula will be useful in the next section. 

6 Topological expansion of the free energy 

We will now compute the topological expansion of the free energy 

oo / rp\ 2h — 2 

^=E(£) *** (6.!) 



h=0 



The free energy of the matrix model satisfy an homogeneity equation. Suppose one scale the parameters of 
our model in the following way 

tf ] -> Ktf ] , T -► kT (6.2) 

the partition function of the model remains unchanged. Indeed the formal integral stays the same due to the 
fact that it only depends on -4^-. Notice that a side effect of this scaling is the corresponding scaling for the 
filling fractions — > Ke j • This means that changing k do not change T 

k^-T = (6.3) 

OK 

Looking at (|6.ip we have that should be homogeneous in k of order 2 — 2h 

= (2 - 2h)T {h) (6.4) 

Equation (|6.4[> can be rewritten as 

J2 = (2 " 2h)J7(h) (6 - 5) 

t a eM a 

6.1 Inversion of the loop insertion operator 

Let us compute here T^ h ' for all h > 1. !F^ ' and will be computed using other techniques. Note that from 
and cu[ h:> are related by 

^--w®*™ (6 - 6) 

so in principle we only have to invert the loop insertion operator. We will propose an operator and act with it 
on uj[ h \ and then check that what we obtain satisfy (|6.6p . 
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6.2 The H- operator 

Following [5J [T2] consider the integral operator acting on a differential (j>(p) (recall the notation of (|4.37[0 : 

t a eM 

As remarked in (|4.38[) the choice 4>{p) = B(p, q) yields — ydx. If we perform an appropriate dissection T> of 
E along the basic cycles and along paths joining the points above infinity then we can write ydx = d'J' for a 
holomorphic "J" on T>. We can then rewrite (|6.7[) by integration by parts (basically a form of Riemann-bilinear 
identity) 

rr=q 

H-4>{-) = - Res 4>{q) / y(r)dx(r) = - Res </>(?)*(?) (6.8) 

q— >Polcs of 4>{q) J r—o *? — ►Poles of <f>(q) 

Again in [T5] the authors proved that 
(2 - ft - 2/>)4% K ) = -4+ibK, •) = Res 4+i (</)*(<?) = E Res (6-9) 

for all (/i, ft) such that 2h + k — 2 > 0, with oii being the branch points (only poles of 4+1 )• Note that the 
h = 0, k = 1 case amounts to (|4.38[) . The proof of this relation is obtained by applying the H- operator into 
the recursive definition of For more details see [T2"] . 

6.3 Computing 

Define the following infinite hierarchy of functions on the moduli space M. 

~ {h) H-JVj.) EjRes^^g) 

2-2h 2-2/i [ ' 

for /i > 1. Our claim is that F^ 1 ' = T^K As a first check we apply dx(p) g y, t to JF^' 1 ', 



(2 - 2^)^^) = -dx(p)^ £ Has *(*)4%) 



(6.11) 



Let us work out the first term 



d f q d 

dx(p) ^(g) = / dx(p) . . t/(r)da;(r) 



= -dS 9 , (p) 
then we have 



dx(p) / 1 



x(p) - x(q) x(p) - x(o) 



i=0 



(6.13) 



-u{ h \p). 
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after moving contours and integrating. The second term of (|6.11[) reads 



by eq. (|6.9|) . Together we have 



dv( P ) 

V Res 9(q)ui h) (q,p) ( 6 - 14 ) 

^ — ' q—HXi 
i 

-(l-2h)4 h \p) 



dx ^W(x)^ h) = -^ h) (6 - 15) 



This is telling us that — J-^ may only depend on T and the filling fractions. 
In fact we can work out any variation with respect to the parameters i a £M, 

( 2 - ^w/ (h) =H*% G£*<p)) ^(p)Eb- *<p) (^(p) 



V Res f f P ^(^(r,-))^^ + V Res *(p)J" (u^fa 



(6.16) 



Ja^) +(l-2/ l )J a (4 /l) ) 

(2-2/i)J a (wf )> 



We are allowed to change the order of the integrations because the integration path never intersect. We use 
(|6.9p to work out the contribution from the second term. 

It is clear now that T^ h ' satisfy homogeneity equation (|6.5p . Indeed, using the previous result we get 

t a eM ° ta t a eM (6.17) 

= -H-ui[ h) {-) = {2-2h)T {h) 

that together with ([6~T5jl . shows that = .T 7 ^, 



3— ^Res^^^Cp) (6.18) 



2-2h^> pem 

i 

with on the branch points and d'I'(p) = y(p)dx(p). 
6.4 Computing JF^ and .T^ 1 ) 

This method determines all terms in the topological expansion of T except for T^' and . These two terms 
must be computed using different techniques. We do so in this section 
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6.4.1 Formula for 

The computation of the planar limit of the free energy is contained in [18j and we report only the final formula 
(only minor notational differences occur). We are seeking a function on the moduli space such that 

- Res ] x k ydx (6.19) 

3 000,002 / 





= (^Res- 

\0 ooi 




= -(fRes 

yd oo 2 




= i> y dx 




= /*«:=(- 




POO! 

= t ydx 

J 00 2 



1 



(6.20) 
(6.21) 



Y i ydx (6.22) 

oo 



(6.23) 



Theorem 6.1 (Theorem 3.2 in |18j) The planar limit of the free energy is given by the formula (refer to 
h4-3ty f or the meaning of the symbols ) 

di d 2 9 2 

2 ^ = E t^ty (y dx ) + E l f J tf (y dx ) + T My ^) + E E e f )j ^ (» dx ) = ( 6 - 24 ) 

k=a j=Q 3 j=l fe=i 3 

= J2 t My (ix ^Y, tatbJa ® Jb ^ ( 6 - 25 ) 

a a,b 

where the notation J a ® Sblp) stands for the integral operator applied to a bi- differential on the two variables. 
The last expression comes using 

6.4.2 Formula for 

The first topological correction to the free energy can be computed adapting already standard methods (see for 
example jTSl [20] ) to our case. We present here the result and derive it in appendix [C] 

^ (1) = -iln^r fl:c 12 nyVi)J (6.26) 

A Loop equations and spectral curve 

In this appendix we derive the master loop equations, that will allow us to write a recursive formula for the 
topological expansion of the Cauchy matrix model. 

One way to find loop equations is by applying infinitesimal changes of variables to the matrix integral, and 
looking at the first order variation of the integral, which must vanish. 

A.l Loop equation: General form 

The loop equations are a glorified version of the following fact: suppose we have a multivariate integral on a 
domain f2, f{x)dx and a vector field x = h(x). Then the action of the vector field on the integrand induces 
the integral of a "total derivative" , namely 
d 



dt 



f(&)d n x = (h ■ V/ + fdw(h)) d n x = / div ( fh) = - / f{x)h{x)d n - x x (A.l) 
n Jn^ ' Jn ^ ' Jon 
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where the last integral is the "flux" of the vector on the boundary of f2. If cither / or h vanish on (or are 
tangential to) the boundary, then the variation is zero and thus yields identities amongst different integrals. In 
our case x are the two matrices and the domain of integration O is the cone of positive matrices. Consider the 
vector field induced by the infinitesimal change of variable 

Mj -> Mj + eSMj . (A.2) 

The partition function Zq is an integral on the space of positive definite matrices; the infinitesimal variation of 
the integrand under (|A.2[) is the divergence of a vector 

Mi -> Mi + eSMi 

(A.3) 

Z c -> Z c + eSZ c + O (e 2 ) , 

In order to have a vanishing variation the vector field must vanish on the set of matrices with at least one 
zero-eigenvalue. Under this condition we will have 8Zc = 0. From the explicit form of Zc (|2.ip . 5Zc = is 
expressed as 

(j(*Mi)) - £ (tr (^MiFi'(Mi))) - (tn (W 2 ( j^f^ ))) = 0. (A.4) 



J (8 Mi) is the factor coming from the Jacobian of the change of variables. This is the trace of the matrix 

9(Afi) fc ,, ' 



For our purposeqj this computation reduces to two simple rules 



J [^tr a \ b\ = tr f A b\ 

\ U-Mi J J U-Mi e-Mi J 

where A and B are in general combinations of constant matrices in Mi and factors of the type ^_ Ml ■ Of course, 
if A and B do contain such resolvent-like terms, the Leibniz rule applies to A and B. 

A.2 Master Loop Equations 

In the Cauchy model two loop equations appear to be crucial. One of them is quadratic in the resolvents and 
the other is cubic. In this section we derive both of them before extracting all information in the following 
sections. 

A. 2.1 Quadratic Master Loop Equation 

Consider the following change of variables for Mi 

SMi = ^=r( — . (A.7) 



8 we leave it to the reader to check these rules. Also, see [8]. 
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Note that the the vector field above is permissible since it preserves the cone of positive matrices, 
any zero eigenvalue of M\ is invariant under the flow; the eigenvalues of M' = M + e^r ( gjjg " 

m- = wij + e|j ( g-m- — f ) >° 0. Using the rules (|A.6|) we find 

J(5M X ) = ( Jtr 



The associated loop equation is 

(wmo) 2 ) - ^/(e) (^i(o) + (mo) + (bi(o) + 7 (ji) = °- 



where we have defined 



T / 1 



s i(0 = -iT7otri -— -tra 



TV 2 1 \Z- M { \M 1 ®I + I®M 2 
I 1 = Itr(^(M 1 )) + ^tr li2 / 



TV v 1V " TV 2 ' \M 1 ®I + I(g)M 2 

Make a similar variation for M 2 

T 2 ( 1 1 

<5M 2 = 



TV 2 \r\~M 2 i] 
with the Jacobian 

( T ( 1 
J(5M 2 ) = ( — tr 

with the associated loop equation 



TV \t)-M 2 



and the definitions 



(W 2 (- V )Y) V 2 \r,) (W 2 {-r,)) + (P 2 (-r?)) + (S 2 (-r/)) + — (/ 2 ) = 



W 2 (-»/) = 5tr / 



TV \rj-M 2 

TV \ r] — M 2 
„ T^_ ( 1 / 1_ 



I t r(^(M 2 )) + ^tr 1 , 2 / 



^TV v ^ v " TV 2 ' V M i® I + I ® M 2 > 
Combining both variations, setting r\ = — £ = m and using the relation 

S 1 (a;) + S 2 (x) = Fy 1 (a;)W 2 (a;) 

we obtain the quadratic master loop equation 

{W^x)) 2 ) + ( (W 2 (x)) 2 ) + (Wi{x)W 2 {x)) - V((x) (Wi(x)) - VZ(-x) <V 2 (x)) 

p 1 (,)) + (p 2 (,)) + i((/ 1 ) + (/ 2 ))=o. 
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Using the functions Yi defined in 13. II the equation (|A.16[) is expressed in a much more compact way 

( (Y 1 (x)f) + ( (Y 2 (x)) 2 ) + (Y 1 (x)Y 2 (x)) = (R(x)) ~\{h + h) ■= (R(x)) (A.17) 

where 

R{x) = \{V{{xf + V 2 \-xf - V[{x)Vi{-x)) - P 1 (x) P 2 (x) 
- {U[(x) f + {U' 2 {x)f + U[(x)U^(x) - P^x) - P 2 (x) 
is, in the limit N — > oo, the same R(x) found in pQ. 

Remark: Note that in [l] the ^ term in the RHS of 1 13, 81 is absent due to the constraints used on the potential which guaranteed 
the support of the equilibrium measures to be disjoint from 0. In this case one can verify that the expectation of I\ + 12 is a total 
derivative vanishing on the boundary of the integration and hence it is zero. 

Remark: Note that R(x) is a rational function in the variable x with only one pole at x = of order 2 coming from the squares 
of Yi. Clearly, although each term on the LHS has cuts along the real axis R, their combination is rational function with only one 
pole at ^ = of order 2 in general, or of order 1 in case t'_{ = 0. 

A. 2. 2 Cubic Master Loop Equation 

For the cubic master loop equation the procedure is exactly the same but with slightly more complicated changes 
of variables. We are going to perform one change of variables for each Mi and M 2 , and then combine them. 
Consider the following change of variables 8M\ and its related Jacobian J (5 Mi) 

iM > - w {jha - \) (h (fte) - h fans) - v! ® + 
asmj = (I, (*))' - (** f^V I* (*) - «'«) - vim (l« 1 



N \£-M 1 JJ \N \t-MxJJ N \ij-M 2 J v iS "'\N \£ - M x 

T 3 ( 1 \ T 3 1 / 1 \ 
N 3 ^ l(C-Mi) 3 J ~ N 3 I (£-Mi) 2 J 



(A.19) 



after a little bit of algebra one gets the following loop equation 

-2 /1 J2 



(2V{(0 viin)) (^i(C) 2 ) + v((0 (wxftWti-v)) 
+ (Si(o (WtiO - w 2 (-v) - v{(0 + vfa))) 

+ ( (pi(0 + jh) (Wi(0 - w 2 (- v ) v{(0 + v^))\ 



(A.20) 



Consider also the symmetric change of variable for M> 



2 



T 3 I 1 \ T 3 1 ( 1 
~N 3 ~ tl \ { V -M 2 f) ~ N 3 ^ I ( V -M 2 f 



(A.21) 
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to obtain 

T 2 ( 1 d 



+ (V{(0 - 2^(77)) (W2(-r?) 2 \ + ^(77) (w 2 (-»?)Wi(0 



(A.22) 

- (s 2 (0 (w 2 (-r?) - - y 2 > ? ) + 

+ ( (AH?) + — /a ) (WaH) - Wi(0 - ^(77) + V\'(£)) 



-»? 

Now, set —77 = £ = x, subtract equation (|A.22|) from (|A.20|) . use (|A.15[) and introduce Yi and Y2 to obtain 
the cubic master loop equation 

(Fx(x) 3 ) - (y 2 (x) 3 ) - (i?(x)(Yi(x) - Y 2 (x))) ^ + 1 A) (^i(x) - W a (s)) = (A.23) 

This equation implies that both Y\ (x) and Y2 (x) satisfy the same equation 

(y fc (x)3) - (*(x)Y fc (x)} - £ (i^ + i A) (W k{x) ) = , for k =l,2 (A.24) 

where -D(x) is still unknown. To find an expression for D(x) we use the same basic variations above to obtain 
the following equation 

(y (x) 3 ) - (r(x)y (x)) ^ QA, + 1A) ( Wo(e) ) = - (y l( x) 2 r 2 (x)) - <Vi(x)y 2 (x) 2 ) (A.25) 

where we define 

W (*) = -Wi(z)-W2(x) 

(A.26) 

F (x) = -Yi[x) - Y 2 (x) = U' [x) - W (x) 
On the other side, from the two cubic master loop equations (|A.24[) we have 

(y (x) 3 ) - (r(x)y (x)) - ^ g A, + 1 A) (m/ (x)) = -2 (i?(x)) - 3 (f!(x) 2 f 2 (x)) - 3 (r!(x)r 2 (x) 2 ) 

(A.27) 

The RHS of equations (|A.25[) and (|A.27|) determine, after some more algebra and more loop equations, 

(d(x)) = - (f!(x) 2 f 2 (x) + yx(x)y 2 (x) 2 ) 

= -U^(x)U{(xM(x) - U{(x) (P 2 (x)} - U' 2 {x) (Pi(x)) - 1 (t/{(x) (l 2 ) + (7 2 (li)) (A.28) 
+ (<9i(x) + Q 2 (x)) + (Wi(x) 2 W 2 (x)) + (^i(z)^ 2 (x) 2 ) + (Si(x)) + (&(*)) 

with 



A 2 V x-Mi V M i« lII + II( 8' M 2 

T 2 (Vi{x)-Vi{M x )^ ( 1 



Q 2 (-x) =--tr — - tr 



A 2 V x-Mi 1 VMi®I + I(8)Af2 // 

T 2 /V?(Af x ) / 1 ( } 



c /j,"! fj. I i v tr I 

U ; A 2 \x-lVh 2 \M 1 ®l + l®M 2 

&H0 =772" tr - — rptri 



A 2 V a; - Mi V Mi <g) I + I ® M 2 
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In [T] the authors already considered the large N limit of this equation and found the large TV limit of (^D(x)^ 
that we call D^°'(x). It is clear that the last 4 terms in (|A.28[) only contribute to the pole structure at x = 0. 
The analysis in [T] shows that the pole coming from these terms is at most of order 2. The actual order of the 
pole at x = depends on t_[ through the products of U[{x). Whenever both t_i > are different, the pole 
is of order 3 and comes from the product Uq(x)U[(x)U2(x) . Whenever the symmetry = is present the 
pole is of order 2 since Uq does not have a log term. 

Remark: Note that now, eq. HA. 2411 is satisfied by k = 0, 1, 2. 

Remark: Note that, by construction, the LHS of equation J A. 241 ) may have cuts only on R+ for k = 1, while for k = 2 the cuts may 
only appear on R_. Indeed, Yi(x) has poles only the spectrum of Mi, x = i, £ R_|_ , while ¥2 0*0 h as poles only on the spectrum 
of —M2, x = —yi G R_. Thus ( D(x) \ is a meromorphic function of x € C everywhere except at x = where it may have poles. 



B Definitions on the algebraic curve 

In this appendix we will give a brief review of several definitions of objects on the algebraic curve. The notions 
are relatively standard and we refer for example to [2TJ [22] . 



B.l The algebraic curve 

Given an algebraic curve E(x,y) = (irreducible polynomial equation of degree d\ in x and c?2 in y), its points 
(x,y) parametrize the points p of a Ricmann surface S. Choosing x € C U {0} (y E C U {0}) as the base 
coordinate there are di (d\) solutions j/i(x) (xi(y)) of the algebraic curve. Those functions define c?2 x-sheets 
(d± j/-sheets) which form a chart on the Riemann surface. Since, for every point p £ E there is a pair (x,y) 
that solve E(x,y) = we can define two meromorphic functions x(p),y(p) on S. For a generic p = p^ £ S 
there arc di — 1 other points pW, i = 1, . . . , cfe — 1 such that x(p^ ) = x(p). Each one of these points lies in 
a different x-sheet. In general, on these points we have y(p^) ^= y(p^) for all pairs The same definitions 
can be made based on y-sheets but we do not use it in this work. 

B.2 The branch points, branch cuts and locally conjugated points 

We will consider curves where dx(p) and dy(p) never vanish at the same point p (i.e. without cusps). The points 
where dx(p) = are called branch points and are denoted by a^. We only consider algebraic curves where the 
zeros of dx(p) arc simple. The sheets are connected pairwise by cuts that go from one branch point to another, 
every branch point having only one attached cut. 

In a neighborhood U a of any branchpoint a there arc two points p,p £ £ such that x{p) = x(p). The map 
: U a 1—* U a is a local involution, namely, p = p. Note that the notion is purely local and has no intrinsic 
meaning at the global level. If one uses the local coordinate z{p) := y/ x{p) — x(a) as a conformal parameter 
for U a then z(p) = — z(p). 

B.3 Cycles and genus 

In every Riemann surface we can choose a basis of cycles for the homology group. These cycles can be chosen to 
form a canonical basis, i.e. two types of cycles Ai,Bi homologically independent which satisfy the intersection 
conditions 

.4, • .4, Bi ■ Bj 

(B.l) 

Ai-Bj = -Bj ■ Ai = Si j. 
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A Ricmann surface of genus g has exactly g pairs of conjugated cycles Ai,Bi. Usually one defines the homotopy 
group by choosing related cycles A and B. These cycles are equivalent to the previous ones with the constraint 
that all pass through a common point po- Cutting £ along A and B defines £ which is isomorphic to an open 
region of the complex plane C. This is called the canonical dissection of £ 

B.4 Abelian differentials. 

There are three type of meromorphic differentials on an algebraic curve. The differentials of the first kind 
(holomorphic) form a vector space of dimension g and we denote them by dut. We can always choose a basis 
dui(p) such that it satisfy 

/ d Ui {p) = 6 itj (B.2) 
With this normalization we define the Ricmann period matrix as the matrix of B periods 

d Uj (p) (B.3) 

Bi 

This matrix is symmetric and its imaginary part is positive definite 3(t) > 0. These differentials dui(p) are 
called normalized Abelian differentials of the first kind. 

The Abelian differentials of the second kind are defined as meromorphic differentials with poles but 
with zero residues. A basis for these differentials may be called dQ^\p) which is a meromorphic differential 
with a pole at p — q without residue normalized as follows 

dn { k q \p) ~ (-kz.ip)-"- 1 + o(i))dz,(p) 
dnp>(p) = o (R4) 

Ai 

in some local parameter z q (p) such that z q (q) = 0. 

Finally the differentials of the third kind arc meromorphic differentials that have only poles of first 
order. Again we will choose a basis for these differentials consisting of differentials d5 g>r (p) with only two poles 
of order one, say at p = q,r, with residue 1,-1 respectively, and with vanishing A cycles, i.e. normalized 
following 

dS,, r (p)~(-l + £)(l))d z ,(p) 

dVb) ~ (—^ + O (1)) dz r (p) (B.5) 



(P) 



j> dS q , r (p)=0 



where the local coordinates are defined as above. The differentials in this basis is called normalized Abelian 
differentials of the third kind. 

B.5 The fundamental bi-differential and relations with the Abelian differentials 

We also introduce the fundamental bi-differential (Bcrgmann kernel) B(j>, q) which is a meromorphic bi- 
differential with poles only at p = q of order two with zero residue, and normalized in the following way 

(*(p)-*00) 2 m 



j B{p,q) 



for i = 1, . . . , g 
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this bi-diffcrcntial is fundamentally connected with the Abelian differentials presented above. 

B(p,q)=B(q,p) 
d/(p) = Res B(p,q)f(q) 

q^p 

dx(q)dx(p) 



*-f (x{p) - x{q)y 

(p B(p,-) = 2iridiij(p) 

dS q Ap) = J B(p, ■) 



(B.7) 



d^\p) = Res B{p,q)z r (q)- k = - Res B(p,q)z r (qY 



q^p 

dS q:r (-) = 2ni / duj(-) 

Bj ' Jr 

t—q /*t—s 

dS s . r {t) = & dS q . p (t) 

I t—p J t—r 

In all these equations, the integrations paths that are not closed cycles are constrained to not cross the A 
and B cycles, in other words belongs to the cut algebraic curve along the cycles. 

C Derivation of J 7 ^ 

For the reader's convenience we derive the first topological correction J^~' following [H?I I20j . 
The computation consists in integrating the relation 

dX ^dk) Til) = ^ (C ' 1} 

The equation (|C.1[) should be taken as a generating function for the derivatives of J-^ with respect to the 
moduli in the potentials with x(p)~ 1 as the formal expansion parameter. The equivalent formula for general 
moduli (say T, 77's and e's) is 

V W = -Ja(foP) (C2) 

with J a read off table ([437]) . 

The differential u>i is obtained from equation 15.41 

The sum extends over all branchpoints, namely, the zeroes of dx(p). The residue around the branchpoint a is 
computed straightforwardly by Taylor-expanding the expression. After a computation in the local parameter 
z a := y/x(p) — x(a) one obtains 



Sb (a) y"'(a) 



96y'(a) ^' ' \24y'{a) 96(y'(a)) 2 / 
W A V- „. > y'"(a) \ (C.4) 



(A) (B) 

We have emphasized the two expression above because each of them is the variation of a separate term, as we 
now explain. We first need the 
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Lemma C.l Let d t denote any variation of the moduli of the curve and recall the correspondence |^.37| ). Let 
a = x{a), dx(a) = be any of the critical values and denote the jet-expansion ofy(z) as y(z) = y(a) + y'(a)z + 
.... Then 

dta = -*J t {B{-,a)) (C.5) 
y'(a) 

= ^J t (fl(^)(,a)-l^fl(,a)) (0.6) 



where the evaluation of the bidifferential and the derivatives are done in the local parameter z = s/x — a. 
Proof. By the discussion preceding (|4.37p we have 

d t y(p)dx(p)=J t (B(.,p)) (C.7) 



Consider the jet expansion of the above identity near p = a in the coordinate z = z a = ^/x — a. Recall that 
the variations are taken at fixed x-value but the critical value a is not fixed and hence dtz = — With this 
in mind we have 

d t y(p)dx(p) =d t (y(z))dx= ({d t y)\ z (z) + y'(z)d t z) dx = 

(d t y)(a) + {d t y'){a)z + (d t y")(a)^ + ••■ - ^ (y'(a) + y"(a)z + ^^z 2 + 2zdz 

( z 1 
= J t lB(;a)+B'(;a)z + B"(;a)— + . 

(C.8) 

where the primes denote derivatives w.r.t. z (applied to the second variable in the case of the bidifferential 
above). Comparing the coefficients of the jet-expansion in z we find 

O t a = -J t (^4) (C.9) 



y'(a) 

^ yik){a) -Jkhy. yik+2){a)dta = TkhY Jt ( Bik+1) ^ a) ) (c - 10) 



(fc + f)r K ' 1 (k + i) 

which immediately yields the assertion. Q.E.D. 

Using Lemma (|C.1I) with k = 1 on the term (B) in (|C.4[) we see that 

w = k dx ^ww ln R y ' {aj) (c - n) 

To identify (A) we need to resort to the results of [55] where (paraphrasing) it was shown that 

Proposition C.l (]22tf) Suppose x : E — > P 1 is a branched covering of the Riemann- sphere with branchpoints 
oij . Then the following differential on the space of critical values of such coverings is closed 

d 1 
d t lnT Bx := -^-(lnT Bx )9taj ■= - ^ S B (aj)d t aj (C.12) 

j 3 

and defines a local function called Bcrgmann tau function of the covering. 

In [23j (and later works by the same authors) explicit expressions for tb x where also derived in terms of 
Thcta-functions. We refer to their work for the explicit expressions since it would lead us too far astray. 
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We can however conclude by identifying the term (A) in (|C.4[) as follows 

d 



(A) = d2;b) ^%) lnT ^ 12 

and hence -finally- 



(C.14) 



One can verify that also the variations with respect to the other moduli (filling fractions, total charge etc.) 
coincide and hence we can set 



(C.15) 
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